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ABSTRACT 

The  propagation  of  plane  waves  in  a  stratified  ionosphere 
is  considered.  It  is  shown  that  the  constmaction  of  'time  •  charac- 
teristic waves  depends  upon  the  diagonalization  of  a  matrix  by  a 
Loewy  transformation.  For  a  homogeneous  ionosphere  this  construc- 
tion can  always  be  carried  out  explicitly.  For  the  general  case 
we  obtain  the  wave  solution  in  terms  of  characteristic  waves  which 
have  the  properties  usually  attributed  to  characteristic  waves  in 
magneto-ionic  theory.  In  addition,  the  interaction  between  these 
waves  is  made  clear  by  the  introduction  of  continuous  reflection 
and  transition  coefficients.  Ordinary  and  extraordinary  wave-re- 
flection levels  are  shown  to  be  points  at  which  the  elementary  di- 
visors of  a  certain  matrix  become  non-simple.  Complicated  expres- 
sions are  obtained  from  which  ionospheric  reflection  and  transmis- 
sion coefficients  may  be  obtained  explicitly. 

Explicit  solutions  are  obtained  for  normally  incident  waves 
\rith   an  oblique  earth's  magnetic  field  and  for  oblique  incident 
waves  with  a  vertical  earth's  magnetic  field. 
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!♦  Introduction 

The  linearized  theory  of  ionospheric  propagation  may  be  separated  into 
three  parts:  determination  of  the  electron  and  ion  distributions,  derivatioi 
of  the  electromagnetic  properties  of  the  ionosphere,  and  solution  of  the  re- 
sulting Maxwell's  equations.  While  all  of  these  problems  have  received  con- 
siderable attention,  the  last  one  has  not  been  advanced  to  a  point  coiqaarable 
to  that  of  the  first  two.  That  is,  the  Maxwell's  equations  which  contain  the 
most  accurately  determined  tensor  properties  of  the  ionosphere  have  not  been 
fully  examined. 

Approximate  solutions  of  Maxwell's  equations  can  be  obtained  by  the  methods 
of  geometrical  optics,  or,  when  separation  of  variables  is  ^plicable,  the  so- 
called  "  wave"  solutions  may  be  determined.  Some  of  the  experimental  ob- 
servations are  not  adequately  accounted  for  by  the  gecmetrical  optics  solu- 
tions. Therefore  we  mast  either  correct  the  geometrical  optics  solutions 
by  including  higher-order  terms  in  the  asymptotic  expansion,  or  seek  accurate 
wave  solutions.  The  latter  alternative  shall  be  considered  in  Part  I  of  the 
present  paper.  The  problem  of  finding  these  wave  solutions  would  entail  first 
an  investigation  of  the  separation  of  Maxwell' s  equations  with  tensors  e,  \i 
and  (T  appropriate  for  the  ionosphere.  However,  we  will  consider  a  horizon- 
tally stratified  medium  and  seek  "plane  wave"  solutions;  for  this  case  the 
problem  can  immediately  be  reduced  to  finding  the  solution  of  a  system  of 
ordinary  differential  equations.  These  separated  equations  can  then  be  solved 
by  the  extended  Peano-Baker  method  described  in  a  companion  report  [2j , 

The  form  of  the  solution  for  an  arbitrary  inhomogeneous  anisotropic 
medium  is  seen  to  justify  the  form  of  solution  which  is  assumed  by  some 
authors  at  the  start  of  their  investigations.  In  the  neighborhood  of  a 
turning  point,  the  solution  indicates  the  extent  to  which  the  four  "  characteristic" 
waves  are  coupled.  From  the  general  formulation  of  the  problem  it  is  easy  to 
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determine  sufficient  conditions  for  the  complete  uncoupling  of  the  » ordinary" 
and  'extraordinary*  waves • 

More  specific  formulations  for  various  combinations  of  the  earth*  s  magnetic 
field  and  incident  wave  directions  are  considered  in  Part  II,  Two  of  these 
cases,  (l)  normal  incidence  with  an  oblique  magnetic  field,  and  (2)  oblique 
incidence  with  a  vertical  magnetic  field,  are  then  solved  by  the  method  of 
analysis  given  in  Part  I, 

PART  I.  GEJIERAL  lOTOSPHERIC  THEOHI 
2»  Formulation 

We  consider  a  plane  ionosphere  with  horizontally  stratified  properties 
in  the  layer  z  <  z  <  z- ,  The  exterior  of  this  region  is  aosumed  to  be  free 
space  extending  to  infinity  in  both  directions.  If  we  suppress  the  time 
factor  e"   ,  Maxwell's  equations  for  this  composite  medium  become 

V  X  E(x,y,z)  =  i  ^  H(x,y,z) 

V  X  H(x,y,z)  =  -i  I  K(z)  E(x,y,z), 

Here  the  magnetic  permeability  has  been  taken  as  unity  throughout.  The  space 
and  the  tensor  dielectric  properties,  K(z),  are  as  follows. 


For  the  ionosphere  (z  <z<z,): 


2  2  ' 

1  -  A(6  -y^)    A(y^y2+i5y^)     ACy^y^-iByg) 


(2.2)  TCz)-  (k^^(z))  a    I   A(y^y2-i5y^)    l-A(5^-y^)      A(y2y3+i5y3_)  |  . 

A(yiy3+i6y2)    ACygy^-iSy^)    l-A(5^-y^) 

For  fJree  space  (z  <  z  and  z  >  z,) 

(2.3)  K(z)  «  I,  the  unit  matrix. 

The  derivc.tion  of  the  above  dielectric  tensor  (2.2)  is  well  known  and  has  been 
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carried  out  in  a  rigorous  form  by  van  der  Wyck  [3} .  The  ionospheric  model 
assumed  is  that  of  free  space  containing  a  distribution,  N(z),  of  free  electrons 
with  charge  e  and  mass  m.  These  electrons  collide  with  heavier  particles  with 
average  collisional  frequency  v(z)  and  are  in  the  presence  of  an  external 
static  magnetic  field  Hp(z)  due  to  the  earth.  The  notation  used  in  equation 
(2.2)  is  related  as  follows  to  the  above  quantities  and  to  the  usual  notation: 
we  have  defined 

5(z)  «  1  -  i^r'S  1  -  i  Z(z) 

2v 


,/  \   Une  N(z)      1     w„\ 
A(z)  s ^  2 — 5~  =  ^^^^ 


m  03-    6{6'^-J'^)  "  ■  (1-iZ)  [(l-iZ)^-I^J 

The  major  problem  is  to  find  the  solution  of  Maxwell's  equations  (2.1) 
which  has  a  specified  behavior  at  infinity  namely,  in  general,  that  it  be  an 
outgoing  wave,  as  is  insured  by  the  Somraerfeld  radiation  condition.  This 
field  arises  from  a  given  source  or  distribution  of  sources  (usually  a  dipole 
specified  by  a  condition  analogous  to  Sommerfeld's).  The  discontinuity  con- 
ditions for  electromagnetic  fields  must  be  satisfied  at  the  planes  z  =  z 
and  z  =  z, .  Such  a  solution  may  be  constructed  by  first  obtaining  plane  wave 
solutions  of  equation  (2,1)  wliich  have  the  required  continuity  and  outgoing 
properties,  and  the  integrating  over  all  possible  orientations  to  produce  the 
required  source.  In  the  present  paper  our  main  problem  will  be  that  of  ob- 
taining such  plane  wave  solutions. 

If  cylindrical  coordinates  are  used  in  the  previous  formulstion,  equations 
(2,1)  and  (2,3)  remain  unchanged  but  the  matrix  (k. .(z))  of  equation  (2.5)  is 
altered  (see  [3],  Chap,  V),  In  this  case  cylindrical  wave  solutions  would 
be  sought  and  they  could  be  superposed  to  construct  the  field  due  to  an  ar- 
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bitrary  source,  A  cursory  examination  of  this  problem  has  shown  that  only 
in  the  case  of  vertical  earth's  magnetic  field  may  cylindrical  wave  solutions 
be  obtained  by  separation  of  variables. 


3»  Generalized  Transmission  Line  Equations 

The  electromagnetic  field  components  which  satisfy  equations  (2.1)  are 
divided  into  longitudinal  and  transverse  fields  by  calling  the  direction  of 
increasing  z  the  "transmission"  direction.  Thus  E-(x,y,z)  and  H^(x,y,z), 
the  z-components  of  the  electromagnetic  field,  comprise  the  longitudinal 
field  and  the  remaining  components  are  the  transverse  field.  From  Maxwell's 
equations  (2,1),  the  longitudinal  field  may  be  obtained  algebraically  in 
tenns  of  the  transverse  field  and  some  transverse  derivatives,  Tliis  yields 
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In  the  remaining  Maxwell  equations  the  longitudinal  field  may  be  eliminated 
by  using  the  above  expressionsj  thus. 
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This  linear  second-order  system  of  partial  differential  equations  involves 
only  the  transverse  field,  and  from  their  solution  the  longitudinal  field 
can  be  computed  using  equations  (3»l)»  Applying  the  method  of  separation 
of  variables  to  the  system  (3.2)  it  is  found  that  there  exist  solutions  of 
the  form 


.  0) 


(3.3) 


-^       --?    i  7(px+qy) 
E(x,y,z)  -  V(2)  e 

_^       ^    i  ^(px+qy) 
H(x,y,z)  -  I(z)  e 


These  solutions  miay  be  interpreted  as  plane  waves  whose  orientation  is  deter- 
mined by  the  arbitrary  parajneters  p  and  q.  Alternatively,  by  a  slight  ex- 
tension of  ordinary  terminology,  we  may  say  that  a  particular  choice  of  p 
and  q  determines  a  particular  "mode"  *  Of  course  in  the  present  case  there 
is  a  continuum  of  modes,  rather  than  a  discrete  set,  as  occur  in  some  elec- 
tromagnetic problens. 

The  equations  for  the  determination  of  V(z)  and  I(z)  are  obtained  by 
substituting  equations  (3.3)  in  equations  (3.1)  and  (3.2): 


(3.U) 


and 


(3.5) 


1^(2)  =   [i  l^{z)   -  qV^(z)] 


.   dio 
ic  2 

CO  dz 
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These  equations  may  be  considered  generalized  transmission  line  equations 
(see  the  treatment  of  Maxwell's  equation  in  isotropic  media  [U]).     As  such, 
they  iBpresent  a  three-wire  non-uniform  transmission  line  whose  distributed 
parameters  depend  upon  the  frequency  of  the  impressed  voltage. 

The  problem  formulated  in  Section  2  has  now  been  reduced  to  that  of 
solving  the  transmission  line  equations  (3»5)  with  the  appropriate  k.  .(z) 
(given  by  equations  (2.2)   and  (2.3))  in  each  region.     The  discontinuity  con- 
ditions associated  with  Maxwell's  equations  are  satisfied  in  the  present 

case  if  the  transverse  field  components  are  continuous,  i.e.,  if 

r-  -1    z=a+o 

(3.6)  L^l(2)»V^^'V^^»V^^J   z-a-o    "  °* 

Here  z  =  a  is  any  level  at  which  there  is  a  discontinuity  of  the  medium  and 
where 

&(x)]^^  £  f(a^)  -  f(a  ). 

The  Sommerfeld  radiation  condition  which  is  to  apply  in  the  free  space 
above  the  ionosphere,   z  >  z, ,  becomes 


av.(z) 

Bz       ~  "  c  "+' j' 


lija     I  -^^ i^XV.(z)|   -  0 

'        Bz  c     +   1        ' 


z->oo 

(3.7)  ,     J  =  1,2. 

BI.(z) 

a->oo  '' 

Here  \  is  the  as  yet  unspecified  propagation  constant.  This  condition  ensures 
that  at  infinity  the  field  components  behave  like 

exp  (  i  ^'^"^  ^^ 
and  thus  represent  upgoing  w;<ves. 

Below  the  ionosphere,  z  <  z  ,  we  shall  not  impose  a  radiation  conditicn 
but,  rather,  specify  the  incident  electric  or  magnetic  field  amplitudes  at 
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CO 


some  level,  say  z  =  0.  Thus  for  j  "  1  and  2  we  require  either 

(3.8a) 

or 

(3.8b) 


av.(2) 


z-O 


=  2i  ?  X.A^  , 


t 


z=»0 


2i  7  X  B,  . 
c  -  J 


Here  X  is  again  the  propagation  constant  and  the  A.  and  B.  are  the  amplitudes 
of  the  transverse  waves  incident  upon  the  ionosphere.  The  conditions  (3.7) 
and  (3.8a)  or  (3.8b)  together  with  the  continuity  requirement  are  sufficient 
to  ensure  uniqueness  of  the  solution.  The  mode  synthesis  could  proceed  only 
after  the  problem  forniilated  above  has  been  solved. 

h»     Solution  of  Transmission  Line  Equations 

The  analysis  is  simplified  by  introducing  the  column  and  square  matrices 
¥(2)  and  A(z)  respectively,  defined  by 
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The  generalized  transmission  line  equations  (3.5)  can  now  be  written  in  matrix 
form  as 
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(U.2) 


dz 


W(z)  -  i^  A(z)W(z). 
c 


This  system  is  solved  by  the  extended  Peano-Baker  method  descrLbed  in  a  pre- 
vious report  [2].  ■&  sumriary  of  the  definitions  and  identities  necessary 
are  included  for  convenience  in  Appendix  I  of  the  present  report.  The  solu- 
tion of  the  system  (U.2)  is 

W(z)  -  il  I    A(x) 


(1;.3) 


J  ° 


in  any  of  the  regions  described  by  equations  (2,2)  and  (2,3).  Here  K   and  z 
must  lie  in  the  same  region  and  W  is  the  value  of  W(g),  In  the  presert 
section  we  shall  confine  our  attention  to  the  ionosphere,  taking  ^  =  z  and 
W  "  W(z  ),  The  solutions  for  other  regions  are  considered  in  Section  6» 
The  wave  behavior  and  other  expected  features  of  the  solution  are  not 
contained  in  the  present  form  (U,3).  To  remedy  this  situation  we  introduce, 
four  new  dependent  variables  u.(z),  defined  as 


(ii.Ii) 


s  U(z)  =  P"'^(z)W(z),     or 


W(z)  «  P(z)U(z), 


.-1/ 


Here  P(z)  is  an  arbitrary  non-singular  hxh  square  matrix  with  inverse  P"  (z) 


and  elements  p.  .(z)j  i,j  =  1,2,3,U,     Thus  the  u.(z)  are  linear  combinations, 
with  non-constant  coefficients,  of  the  w.Cz),   and  vice  versa.     From  equations 
(U.3)   and  (U.lj)  the  equation  satisfied  by  U(z)  is  derived: 


(U,5) 


^  =   [i  2  p-^(z)A(z)P(z)-  P-^(Z)  ^]  U(z), 


This  system  has  the  unique  solution 


(ii.6)        U(z)  =  -H-  - 
where 

However,  since  P(2)  is  arbitrary  we  may  attempt  to  choose  it  such  that  the 
solution  ih*6)   has  some  of  the  desired  features. 

If  two  square  matrices  L(z)  and  M(z)  are  related  such  that  for  some 
non-singular  matrix  P(z), 

P'^(z)L(z)P(z)  -  P"^(2)f(z)  -  M(z), 

then  they  are  said  to  be  equivalent  in  the  sense  of  Loewy  [6],  Thus  if  i  —  A(z) 
is  equivalent  in  the  sense  of  Loewy  to  a  diagonal  matrix  J\.{z)   then  from  (li.6) 
and  (1.2)  (see  Appendix  I)  we  see  that  each  u.(z)  would  be  the  exponential 
integral  of  the  i-th  diagonal  element  of  J\{z),     Such  a  solution  represents 
true  characteristic  waves  whose  existence  has  been  considered  by  Booker  [$] 
and  others.  If  the  medium  were  homogeneous,  A  and  P  woaild  be  constant  and 
hence  P  would  equal  Oj  thus  the  Loewy  transformation  would  reduce  to  the 
ordinary  similarity  transformation  P~  AP,  In  this  case  the  conditions  for 
diagonalization  are  well  known  [7]  and  indeed  diagonal ization  may  be  carried 
out  for  various  ionospheric  configuration.  This  is  in  conformity  with  the 
accepted  theory  that  'true* characteristic  waves  may  exist  in  a  homogeneous 
ionosphere.  However,  in  the  general  ionospheric  case  the  existence  of  the 
diagonalizing  Loewy  transformation  and  hence  the  existence  of  true'  charac- 
teristic waves  remain  open  questions. 

Perhaps  the  next  best  simplifying  choice  of  P(z),  as  suggested  above 
for  the  homogeneous  case,  is  that  in  which 

(U.7)         p-^(2)A(z)P(z)  -  A(z)  z   (>i(z)6^p, 

which  is  a  diagonal  matrix.  Such  a  P(z)  always  exists  and  can  be  easily 
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obtained  if  the  eigenvalues,  X.(z),  of  A(z)  are  distinct,  i.e.,  if  there  are 
no  repeated  roots  of  the  characteristic  equation 

(1|.8)  |U  -  A(z)|  =  0. 

This  algebraic  equation  will  be  shown  to  be  the  well-known  fourth-degree 
equation  of  magneto-ionic  theory.     The  inequality  of  the  roots  of  this  equa- 
tion is  a  sufficient  but  not  necessary  condition  for  (li.7)  to  hold.     The  more 
general  requirement  is  that  the  matrix  (\I-A(z))  have  simple  elementary  di- 
visors  [i]  •     This  distinction  will  be  recalled  in  Section  6  where  it  is  of 
importancej  but  for  the  present  we  assume  the  eigenvalues  to  be  distinct. 
The  solutions  u.(z)  given  by  (U.6)  and  corresponding  to  P(z)  takem  to 
satisfy  (i|.7),  where  the  eigenvalues  are  distinct,  are  called  characteristic 
waves. 

To  examine  the  fonii  of  these  characteristic  waves  we  introduce  the 
notation 

-1/ 


(U.9)  (^ii^z))  5  -  F     (z)  P(z). 


Next,  the  coefficient  matrix  of  the  system  (U.5)  is  split  into  a  diagonal 
matrix  and  a  zero-diagonal  matrix  as  follows: 

(il.lO)  D(z)+N(z)  s   [i>"^(z)A(z)P(2)-F"^(z)f(z)l, 

where  the  diagonal  matrix  is 

(li.lla)       D(z)  1  (li  I  X>r^^]6.j^) 

and  the  zero-diagonal  matrix  is 

(it.llb)  N(z)  z   iT.^iz)\l-b.^). 

With  the  above  definitions  the  solution  (Ii.6)  becomes 
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U(z)  -  -0-    D+N  1  U 


(1|.12) 


exp 


.].[ 


y   D(x)dx  H  <  expj-/ D(€)di^  N(x)  exp 


/D(€)d? 


^^o 


U 


(see  Appendix  I  for  the  identities  used).  To  simplify  further,  we  intro- 
duce the  matrix  M(x)  with  elements 


(ii 


f 

.13)     ni^.(x)  -  [1-6^ .]  |r^ .(x)exp 


t 


/  [r,,U)-r ,,(«)]<« 


-^ 


'  exp 


i.'^/[\U)-\^U)]dK 


This  matrix  is  the  argument  of  the  matrizant  of  the  final  equation  (Ii.l2)  which 

has  been  written  explicitly  by  means  of  equations   (it. 11)   and  (I»2)  and  (I.I|)« 

Letting  the  elements  of  this  matrizant  be    -*!..  .(z)  we  have,   from  equations 

(i;.12),    (U.13)   and  (I.l), 

z  U      z  ^2 

(U.IU)       -^ij(z)  2  ^ij+Z'^^^ij^^^*  11/^2/ ^"^ik^^2Hj^^^'^   ***• 

\  k=I  z^         Zq 

Finally  we  introduce  the  column  matrix  C(z  ,z)  with  the  four  elements  c.(z  ,z) 

defined  by 

h 
c^(z^,z)  zTl  -n-i;5(z)^j(Zo^*  ^  "  1*2,3,1*. 


(Ii.l5) 


From  equations  (li.l2)   to   (U.l5)  we  obtain  for  the  individual  characteristic 
waves 

(1;.16)  u.(z)  =  exp(/r^^(x)dx).  expQi  ^  A   (x)d^.  c  (z^,z)j 


3  -  l,2,3,Ii. 


ksi.     Interpretation  of  Characteristic  Wave  Solution,  Each  characteristic 
wave  (ii.l6)  consists  of  three  factors.  The  second  factors  in  those  four 
waves  are 


exp[i  ^/  X-(x)dxJ:    j  =  l,2,3,itj 
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these  reveal  the  *wavelike'  behavior  of  the  characteristic  waves  and  thus  are 
callod  the  "wave  factors*.  They  are  all  distinct  since  the  eigenvalues  X.(x) 
have  been  assumed  unequal.  Let  two  of  these  factors,  say  those  containir^ 
X- (x)  and  \,(x)  represent  rising  waves  and  the  remaining  two  downcoming  waves. 
The  X^(x)  and  X^Cx)  wave  factors  shall  be  called  'ordinary  wave  factors" 
and  the  remaining  two  "extraordinary  wave  factors".   Thus  u,(z)  ani  UpCz) 
are  the  rising  and  the  downcoming  ordinary  characteristic  waves  respectively, 
while  u,(z)  and  111(2)  are  the  rising  the  downcoming  extraordinary  character- 
istic waves  respectively.  This  is  the  usual  ionospheric  terminology  which 
is  justified  in  Sections  9  and  10,  where  specific  cases  are  treated. 
The  first  factor  of  the  solution  CU.I6), 


I  /  r,,(x)dx|,  i   =l,2,3,ii. 


exp,  /  4.  . . 

Jz  ^^ 
0 

yields  a  W.K.B.  behavior  of  the  characteristic  waves.  This  is  shown  explicit- 
ly for  the  troposphere  in  [2]  and  for  the  two  specific  ionospheric  cases  in 

Sections  9  and  10.  In  all  three  cases  it  is  seen  that  r..(x)  contains  the 

1/2 
negative  logarithmic  derivative  of  X/  (x).  In  the  more  general  case  it  has 

been  shown  (see  |_2j ,  Appendix  III)  that  some  such  logarithmic  derivative  is 

always  present  and  thus  represents  a  generalized  W.K.B,  behavior. 

The  factors  c.(z  ,z)  in  the  solution  are  called  'continuous  scattering 

factors" .  They  describe  the  continuous  multiple  reflection  which  is  ejqjected 

in  an  inhomogeneous  medium.  Besides  this  exi^ected  phenomenon,  the  anisotropy 

causes  a  change  of  type  (i.e.,  ordinary  into  extraordinary  and  vice  versa) 

oflection  as  well  as  during  continuous  transmission.  This  change  of 

lich  occurs  in  the  direction  of  propagation  shall  be  called  a  ''  transi- 

tic      iquations  (1;.13)  and  (U.lU)  indicate  that  ll_..(z)  represents  the 

total  contribution  to  a  X.-wave  factor,  by  a  X^-wave  factor,  by  all  possible 

J-  %j 
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multiple  continuous  reflections  and  transitions  and  combinations  of  the  two. 
More  explicitly  in  equation  (U.13)  we  may  introduce 

(U.17)       Rij(x)  •  r.j(x)  exp[/ir^.U)'r..(K))d{\,  ±  fl  3, 

which  are  interpreted  as  continuous  reflection  and  transition  coefficients 
for  X. -waves  changing  into  \. -waves.  From  the  previous  description  of  the 
characteristic  wave  factors,  the  above  coefficients  may  be  classified  as 
follows : 

Inflection  Coefficients     Transition  Coefficients 

No  Type  Change:    ^2*^21*  ^3U»\3 
(U.18) 

Change  of  I^pe:    Rii^jR^J  ^23*^32     ^3'^*  ^2U'\2 

This  table  together  with  Equation  (U.l?)  constitutes  an  extension  to  aniso- 
tropic media  of  the  usual  Fresnel  coefficients  for  isotropic  media.  Such 
an  extension  has  been  mentioned  by  Breraner  |8J,  but  who  does  not,  however, 
give  explicit  expressionsj  this  is  done  here  in  Sections  9  and  10  for  two 
specific  ionospheric  cases. 

As  an  example  of  the  above  description  consider,  say,  J'L.,(z),  the 
contributions  to  the  downcoming  extraordinary  wave,  (\i -wave)  originating 
from  the  rising  ordinary  wave  (\^-wave).  From  Equation  (ii.lU)  one  of  the 
fourth-order  terras  in  this  sum  is 

_  3Ci  X^  Xp 

.-0_^(z)=  ...  +/dx^/  dx    /   dxg/   dx^m^3(x^)m3^(x3)m^2^X2)m2i(x^)+   • 

z  z  z  z 

0000 

( 't  Fourth- order™  means  that  the  number  of  continuous  reflection  and  transitions 
totals  four.)  From  Equations  (U.13)  and  (U.l?)  this  tern  is  observed  to  con- 
tain the  reflection  and  transition  coefficients  Ri  ,,  Hot  »  R-i  o  ^^^   ^21  '^"  that 
order.  This  accounts  for  the  follovring  phenomena  (see  Figure  l):  the  rising 
ordinary  wave  feeds  by  reflection  the  downgoing  ordinary  wavej  this  in  turn 
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feeds  the  rising  ordinary  waves;  and  ty  traasition  this  feeds  the  rising  extrar- 
ordinary  wave,  which  finally  feeds,  "by  reflection,  the  downgoing  extraordinary  wave. 


X.-  wave 
(ordinary) 

AAAAAA* 


X,-wave 


Xj-wave 


(ordinary)    (extraordinary) 


AAAAA 


R 


21 


VV¥V»» 


AAAA 


^43 


^31 


R 


12 


Xg'wave 
(ordinary) 


X4-wave 
(extraordinary) 


Figure  1, 


The  ordinary  and  extraordinary  waves  are  said  to  be  coupled  if  a  wave 
of  one  type  feeds  or  produces  a  wave  of  the  other  tjrpe,  that  is,  if  an  in- 
dependently excited  ordinary  wave  would  give  rise  to  an  extraordinary  wave 
ani  vice  versa.  In  general,  this  is  the  case  for  the  ionosphere.  However, 
if  those  continuous  reflection  and  transition  coefficients  which  refer  to 
a  change  of  type  were  zero,  there  would  be  complete  uncoupling*  From 
Equation  (U.l?)  and  the  table  (U.18),  the  conditions  for  complete  uncoupling 
are  found  to  be 
(U.19)       r^i^  -  r^  -  r23  =  r^^  °  ^13  '  ""31  "   ""sU  "   ''U2  "  °- 

More  generally  it  is  seen  from  Equations  (U.6)  and  (li.9)  that  if  P(z)  is 
chosen  such  that 
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(U.20)        (^  P-^AP-P'^)- 


-   0   0 
-00 

0   0- 

0   0- 


where  the  dashes  represent  arbitrary  functions,  then  there  is  uncoupling. 
In  this  case  the  fourth-order  system  reduces  to  two-second-order  systems. 
There  are  some  known  ionospheric  cases  in  which  this  happens  [9] ,  "the  homo- 
geneous ionosphere  being  a  trivial  example. 

Ub.  Polarization  in  the  Ionosphere.  Performing  the  indicated  matrix  mul- 
tiplication in  Equation  (U.U)  gives 

(U.21)       w.(z)  -  >   p.  .(z)u.(z),  i  =  1,2,3,U. 

j=l  ^    ^ 

From  definition  (U.l)  the  w.(z)  are  seen  to  be  propertional  to  the  electro- 
magnetic field  components^  and  hence  these  components  are  linear  combina- 
tions, with  variable  coefficients,  of  the  four  characteristic  waves.  If 
each  characteristic  wave  is  considered  to  belong  to  a  separate  (i.e., 
independently  propagating)  electromagnetic  field  we  may  then  determine  the 
polarization  of  the  characteristic  fields.  Since  w-(z)  and  w.  (z)  are  the 
X  and  y  components  respectively  of  the  electric  field  vector,  and  w,(z)  and 
Wi  (z)  are  the  corresponding  magnetic  field  components,  the  polarizations 

become: 

Ordinary  Wave  Extraordinary  Wave 

(U.22)  Rising       DoT-mcoming       Rising       Downcoming 

B-polarization     V-^iz)/p^{z)       P3^2^z)/p^2(2)    Pl3(z)/Pi^3(z)  Pil^(z)/P]^l^(z) 
H-polarization    P2^(2)/p^^(z)  p^^iz) /p^^(z)  P23(2)P33(z)   P2]^(z)/p^^(z) 
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However,  these  ratios  give  only  the  projections  of  the  polarizations  in  the 
transverse  plane,  since,  in  general,  the  longitudinal  field  components  do 
not  vanish.  Further,  the  assumption  of  independently  propagating  character- 
istic waves  holds  only  for  true  characteristic  waves  and  must  thus  be  re- 
garded as  an  approximation  in  general.  These  "  polarizations"  apply  only 
to  fields  within  the  ionosphere,  and  there  is  no  justification,  at  present, 
for  assuming  that  a  wave  which  has  come  from  the  ionosphere  will  have  a 
polarization  which  is  the  limit  of  one  of  the  above  expressions.  This 
question  is  discussed  further  in  Section  6. 

Although  eight  seemingly  different  polarizations  have  been  defined, 
it  is  expected  that  they  should  reduce  to  at  most  four,  with  two  pairs 
differing  only  in  sign,  ^lis  is  found  to  be  the  case  for  the  problems 
considered  in  Sections  9  and  10,  although  it  has  not  been  proved  in  general. 

5.  Reflection  Levels 

The  solution  of  Section  k   has  been  obtained  under  the  assumption  that 
A(z)  could  be  diagonalized  at  all  points  in  the  interval  z_  <  z  <  z,.  We 
call  a  point  at  which  diagonalization  is  impossible  a  •'reflection  level" 
and  proceed  to  justify  the  terminology  and  examine  the  nature  of  the  solution 
at  such  a  point. 

If  A(z)  cannot  be  diagonalized  then  XI-A(z)  must  have  non-simple  elementary 
divisors,  and  hence  the  characteristic  equation  (U.8)  must  have  repeated  roots. 
Specifically,  let  us  assujne  that  \^(z)  »  X^iz) ,   and  that  the  other  pair 
is  distinct  in  some  small  interval,  say  z'  <  z  <  "  «  It  is  then  possible, 
in  general,  to  find  a  metrix  P(z)  such  that  [?] 
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(5.1)       P'^(z)A(z)P(z) 


0 

0 
0 


X^iz)     0         0 
0       X  (z)  0 


0         0       X^^(z) 

The  procedure  of  Section  ky  Equations  (U.9)  to   (li.l5)j   is  now  repeated  with 
the  above  choice  for  P(z).     The  forraalism  is  identical,  and  the  main  difference 
occurs  in  equations   (U.H)  which  now  become 


(5.2) 


D(z-)  • 


^(z)     0 

0 

0 

0         ^^(z) 

0 

0 

0           0 

\^{z) 

0 

0           0 

0 

x^(.) 

N(2) 


0        T  *^12        ^13 


^31 


^23 


^32 

%2 


^lli 
'3U 


1*3 


/ 


where,  as  before,   (r.  .(z))  »  -P     (z)P(z),     The  m.  .(x)  defined  by  Equation 
(I1.I3)  differ  ejqDlicitly  only  in  nuAx)   and  m^Ax) ,  which  become 


(5.3) 


m, 


21' 


(x)  -  \  r     (x)  exp(  /    [t  JK)'r^Ax)}dK 
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¥e  obtain  J\.  .(z)  and  c.(z  ,z)  by  using  the  new  quantities  in  equations 
(U.lU)  and  (U.l5)  and  replacing  z  by  z'  whenever  it  occurs.  The  ordinary 
characteristic  wave  solutions  are  obtained  by  using  the  matrizant  identities 
and  carrying  out  the  matrix  multiplication  in  equations  (i;.12)  to  (U.l5)  with 
the  changes  noted  above  and  in  equations  (5.2)  and  (5.3).  ^'^s  jd-elds 


u,(z)  =  expCyr^(x)dx)  •  expD- |y*  X^(x)dx)«  c^(z',z) 
z'  z' 

U2(z)  -  expC/r22(x)dx}  •  expfi  ^  /  X^(x)dx)»  c^{z^  ,z). 


The  expressions  for  u_(z)  and  u.  (z)  are  similar  but  contain  distinct  character- 
istic wave  factors. 

Since  the  two  ordinary  wave  factors  of  equations  (5»il)  are  identical,  no 
distinction  can  be  made  in  this  case  as  to  the  rising  or  downcoming  character 
of  the  waves.  It  may  be  said  that  the  two  ordinary  waves  have  become  so 
strongly  coupled  as  to  be  indistinguishable.  Even  if  the  medium  were  homogeneous 
(i.e.,  if  f(z)  »  0  and  hence  r.  .(z)  •  0)  there  irould  be  coupling,  as  may  be  seen 
from  the  N  matrix  of  equation  (5.2).  This  coupling  would  not  involve  the  extra- 
ordinary characteristic  waves. 

In  inhomogeneous  media,  however,  there  is  coupling  of  the  extraordinary 
to  the  ordinary  waves  at  a  reflection  level.  This  seems  to  contradict  some 
of  the  literature  (e,g,,[.10J)  which  states  that  extraordinary  waves  propagate 
unchanged  in  the  neighborhood  of  an  ordinary  wave  reflection  level.  Of  course 
the  coupling  will  be  small  for  slowly  varying  media.  Such  coupling  is  also 
present  in  the  neighborhood  of  an  extraordinary  wave  jreflection  level. 

The  solutions  indicated  in  the  present  section  may  possibly  be  used  to 
connect  the  solutions  above  and  below  the  reflection  region.  An  alternative 
connection  could  be  i\imished  by  the  solution  expressed  in  equation  (U.3)  for 
the  field  components] we  would  then  use  the  identity  (1.5).  A  final  possibility 
is  to  use  the  solution  given  by  equation  (li.l6)  up  to  a  point  near  the  reflection 
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level  and  then  to  take  the  integrals  along  a  path  in  the  complex  z«-plane 
which  avoids  the  reflection  points.  It  should  be  remembered  that  the  above 
connections  refer  to  the  dxact  wave  solutions  and  not  to  their  asjinptotic 
forms.  It  is  the  latter  problem  which  has  received  much  attention  in  the 
literature  [9]  • 

6.  Free  Space  Solution 

Although  the  free-space  plane  wave  solutions  of  Maxwell's  equations 
are  well  known  and  easily  obtained  it  is  necessar;/^  for  our  purposes  to  find 
the  form  they  take  when  the  notation  of  the  previous  Sections  is  usedj  these 
solutions  are  then  used  in  Section  7«  Thus  we  consider  Maxwell's  equations 

(2.1)  with  K(z)  given  by  (2.3).  Letting  the  field  vectors  have  the  foim 
(3«3)>  we  obtain  by  the  analysis  of  Section  3 

7^(2)  -  -p  12(2)  +  ql]^(z) 

(6.1) 

I^(z)  -  pV2(z)  -  qV^(z), 

and 

(6.2)  ^  =i^AW(z), 

where  W(z)  is  defined  in  equations  (ii.l)  and  where 


(6.3)    A  5 


0 

I-P*- 

pq 

0 

.-,^ 

0 

0 

pq 

-pq 

0 

0 

.i.p^ 

0 

-pq 

-l.q^ 

0 

This  matrix  has  the  characteristic  equation 
(6.1i)        j\I-A|  s  [x^-(l-pW)]^  -  0. 
Thus  the  eigenvalues  are 


/T-p" 


X  =  /  1-p^-q' 


^^'^^  \  -^3  --^2  "h  '^J 
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where  the  subscripts  are  chosen  to  conform  with  the  descriptions  given  in 
previous  Sections.     There  are  two  pairs  of  equal  eigenvalues  but  it  may  be 

shown  that  the  elementary  divisors  of  \I-A  are  simple,  provided  that 

2   2 

p  +q  f   1,  We  assume  this  to  be  the  case.  Then  a  diagonalizing  matrix 


F  and  its  inverse  P"  are 


P  2 


(6.6) 


..-1 


1 
2 


/ 


-pq 

\    ° 

-1 


(1-P  ) 


1-P 

-X 

0 
-pq 


-pq\-^(l-p2)-^ 


pq 

0 
X 

-l.q^ 


M   2x-l 

(1-q  ) 


^ 


0 

pq 

-i+p 


2 


,-1,,   2v-l 
pqX  (1-q  ; 


-X 


-1 


pqX"  (l-q*^) 


0 


2x-l 


,-1,^   2x-l 
-pqX  (1-P  ; 

,-1 


(■^       2^-l 
-(1-q  ) 


n   2x-l 

-(1-P  ) 


-X 


-1 


The  solution  of  the  system  (6,2)  is  now  easily  obtained  by  the  extended 
Peano-Baker  method  of  Section  U. 

We  introduce  the  new  variables,  u.(z),  by  the  column  matrix 

(6.7)  U(z)  =  P"-'v(z). 

Since  f  =  0,  this  matrix  satisfies  the  equation 

(6.8)  dUU)  .i^^U(z), 

where  by  equations  (6.3)  and  (6.6)  F'-'^AP  -/Va  (X^6^j).  The 
solution  of  the  system  (6,8)  is 


U(z) 


z 
z„ 


('M 
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U(z^) 


(6.9) 


e         u^(z^) 


-i  f  X(z-z^) 
e  UgCz^) 

e  U3(z^) 

-i  ^  X(z-z^) 


where  we  have  used  the  definitions  (6.5),  and  U(z  )  ■  P""vir(z^)j  z^  is  an 
arbitrary  point  at  which  the  field  is  assumed  to  be  known.  The  field  quantities 
W(z)  are,  from  equations  (6,7)  and  (6,9), 


W(z)  =  PU(z) 


(6.10) 


U(z^,). 


This  is  the  usual  plane  wave  solution  for  free   space. 

The  coefficient  matrix  (6,3),  indicates  that  if  p  =  0  or  q  =  0   (i,e.       if 
the  wave  is  normal  in  the  xy-  or  the  xz-  plane),  then  V,(z)  and  l2(z)  are  in- 
dependent of  Vp(z)   and  I, (z).     These  two  pairs  of  uncoupled  field  quantities 
are  the   familiar  E-waves  and  H-waves,  which  can  be  excited  independently.     From 
the  discussion  of  coupling  in  Section  U,  it  is  seen  that  the  ordinary  and  extra- 
ordinary waves  are  the   analogues,  in  the  ionosphere,  of  the  usual  E  and  H  waves 
in  isotropic  media. 


7.     Ionospheric  Reflection  Coefficients 

In  order  to  find  the  ionospheric  reflection  coefficients  we  must  solve  the 
problem  formulated  after  equations   (2,U),     To  suirjnarize  briefly,  equations  (2,1) 
and   (2.3)   are  to  be  satisfied  for  z  >  z,    and  z  <  z     >  0  and  equations  (2,1)   and 
(2.2)   are  to  be   satisfied  in  the  range  z     <  z  <  z-,.     By  separation  of  variables 
and  the  reductions  indicated  in  equations  (3.1)   to   (3.5)   and  (6.1)   to  (6.3)   this 
problem  is  reduced  to  the  systems  (U,2)  in  the  ionosphere,  z    <  z  <  z-^,  and  (6,2) 
in  free  space,   z  <  z     and  z  >  z, ,     The  continuity  conditions  (3«6)  must  hold 
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for  z  =  z  and  z  =  z, ,  and  the  radiation  and  incident  plane  wave  conditions 

(3.7)  and  (3.8)  respectively  must  hold  in  the  free  space  regions  indicated. 

With  the  above  formulation  we  seek  the  electromagnetic  field  below  the 

ionosphere,  z  <  z  -  and  in  particular  the  amplitudes  of  the  downcoming 

waves  at  z  =  0. 

By  the  analysis  of  Sections  k  and  6  the  column  matrices  W(z)  which 
are  proportional  to  the  transverse  electromagnetic  field  quantities,   are 


found  to  be 


(7.1) 


where 


W_(z)  =  P_U_(z), 
W(z)  =  P(z)U(z), 
W_^(z)   =  P^U_^(z), 


z  <  z 
—     o 


Z      <    Z    <    Zt 
o   —        —      1 


z  >  z. 


U    (z)   =     Jl    i 


u. 


(7.2) 


U(z)   =     XL 


{ 


U    (0) 


1  c 


"A(2)   -  P'^(z)P(z)fu(z   ) 


U^(z) 


z 

-9-  I 


c  ■**+ 


ru^(zi). 


Here   the  -  and  +  subsciripts  refer  to  quantities  below  and  above  the  iono- 
sphere respectively,   and  the  matrices  P(z),^(z),     P  ,   and  J\     are  de- 
fined in  Sections  U  and  6.     Since  both  regions  exterior  to  the  ionosphere 
are  free  space  we  have  P     =  P     and  J\     =  J[^, 

Now  we  apply  the  first  of  the  conditions   (3-6)   to  W  (z)  and  use  equations 
(6.6)   and  (6.9)   (for  z     "0),  and  the  first  of  (7.2)  to  obtain 


■K- 


(7.3) 


Thus  we  have 


X_u^(0)   +  pqu^(O)   =  A^ 
-  pqu^(O)   =  A^ 

/VVx    \ 


(7.U) 


U   (0) 


U2(0) 

-Ag/pq      / 

and  the  quantities  mAS^)   and  u,  (O)   remain  to  be  determined.     The  radiation 
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condition  (3.7)  applied  to  W  (z)  yields,  in  similar  fashion 

■u,(z,)  \ 

0 


(7.5) 


u^(zi) 


U^(2j^) 

0    / 


where  u,(z-)  and  u-,(z,)  are  to  be  determined. 

To  obtain  these  four  \inknown  constants  we  s?)ply  the  continuity  condi- 
tions (3.6)  at  z  =  z     and  z  =  z,: 


(7.6) 


W_(Zq)   =  W(z^) 


W(z^)   =  W_^(z^). 


Using  (7.1)  and  (7.2)  we  obtain  from  the  above  equations; 


P+U^(2-l)  =  W_^(z^)  =  P(z^)u(z^)  -  W(z^) 


=  P(z^)  SI   < 
o 

^1 
=  P(zi)  SI   . 


-1^  \^^ 


i  7A(x)-P"-'(x)P(x) 


h  P"\z^)W(z^) 


iM-P-^  P-^(Z  )P.  i{  i  1^ 


?f- 


U_(0)^ 


or  finally 


(7.7) 


U3(Z3_) 


z  r 

=  rp;^p(z^)  Jlji  ^  A(x)-p"^(x)p(x)  >  p"^(z^)p_ 

0   /  /A^+A^A,^ 

exp(i|yl_(z^))  (  U2(0) 


This  is  a  system  of  linear  algebraic  equations  for  the  determination  of  the 
reflection  and  transmission  coefficients.  The  large  matrix  of  coefficients 
is  a  special  case  of  the  general  connection  matrix  defined  in  [2j,  Section 
IV,  The  system  (7.7)  yields  for  the  reflection  coefficients 
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u,(0)  .  ^^21V^2UV)  4^  -^Y\U-^2U\3>  r,    e^[i2  ^  X  (z   )], 
2  ^^22^UU  ■  ^2UV^  ''      °     -     °  '' 

(7.8)  (A^+A2)  A^ 

..    ,.,   .   ^V^22-^2lV^  -Xl ^\3^33-^23V^  i^  r-o'-W      ^1 

u,  (O;  =  TT- — 7- i;- — r — c *-*-  exp|i2  -  X  (z   )), 

Here  the  b.  .  are  elements  of  the  coefficient  matrix  of  equation  (7,7),     The 
transmission  coefficients  follow  immediately  if  we  substitute  equations  (7,8) 
in  (7.7). 
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Part  II »      SPECIFIC  CASES 

8,     Some  Specific  Formulations 

Electromagnetic  plane  waves  in  the  ionosphere  have  been  shown ^  in 
Sections  2  to  U,  to  be  of  the  form 

^(x,y,z,t)   =  V(z)  exp[-i  ^(ct-px-qy)! 
-♦  c  — 

^(x,y,z,t)   =  l(z)  exp[-i  ^(ct-px-qy)]      . 

The  six  quantities  v(z)   and  l(z)   are  determined  by 

V3(z)   =  5^^^  [pl^{2.)-ql^izhk^^{z)Y^{z)+k^^{z)^^{z)2 
(6.2)  3^ 


I    (z)   =    [p  V^Cz)   -  q  V^(z)] 


and 


(8.3)  Mil    =i^A(z)w(z). 

Q.4J  Q 


The  matrices  W(z),  a(z)   and  the  elements  k.  .(z)   are  defined  in  equations 
(U.l)  and  (2.2)    (with   (2.li))  respectively.     A  general  form  of  the  solution 
of  the  system  (8.3)  has  been  obtained  in  Section  U.     In  the  present  Section 
we  consider    "  explicit  formulations"    of  the  problem  of  obtaining  such  plane 
wave  solutions.     We  consider  two  directions  of  incidence  of  the  plane  wave, 
and  with  each  of  these,    several  orientiations  of  the  earth's  magnetic  field. 
For  each  case,  we  determine  A(z)   as  x^ll  as  its  characteristic  equation  and 
eigenvalues.     The  method  of  Section  U  could  then  be  applied. 

It   should  be  noted  that  the  characteristic  equation  in  each  of  the 
folloijing  cases  is   just  the  usual  fourth-degree  equation  of  magneto-ionic 
theory  for  that  case.     The  squares  of  the  eigenvalues  are  the  so-called 
ordinary  and  extraordinary  dielectric  constants.     The  elements  k^,(z), 
k-p(z)   and  k^-(z),  which  are  required  for  a  complete  determination  of  the 
field,   arS;  not  given  here  but  may  be  obtained  from  equations   (2,2)   and  (2.U). 

In  bi^ackets,   after  the  name  of  each  case,  are  references  to  some  papers 

in  which  this  particular  case  has  been  considered. 

/ 
A.     Plane  Waves  Parallel  to  the  y-axis;    q  =  0 

In  all   such  cases  A(z)  has  the  general  form 
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(8.U)        A(z)   5 


/-a 

M-' 

0 

-c-«- 

a 

-a-"<- 

0 

^;- 

Y 

c 

0 

P' 

0 

0 

-1 

0 

with  the  characteristic  equation 

\XL-A\   =  X^+[a+a-]x^+[aa%-aix'+;3'JX^+[p'(a+a-;0-(YC*+Y^c)]\ 

^  '^^  +  [(p'aaK--acc^0-(YC-"-a*+r^ca)-n')at3'+YY^0]  =  0, 

This  equation  becomes  biquadratic,  and  may  be  easily  solved,  if 

(a.6)  ^  "  ^"-  =  ° 

YC-"<-  +  Y^i-c  =  0 . 

The  eigenvalues  in  this  case  are  given  by 

X^  =  i  (a^+ati'-|3')   +  i  V(a^+aii'-p')^-U[2aYC-"--p'(a^+aii')-(acc-x-+[j,'Yr*)] 
(8.7)  I 

The  elements  and  the  eigenvalues  of  the  matrix  (8.U)  are  given  below  for 
the  indicated  magnetic  field  directions. 

i)  Vertical  Magnetic  Field:  y,  =  y.  =  0       [llj 


1/2 


(8.6) 


a  =  a-x-  =  c  =  C-K-  =  0 
a  =  13  =  1  -  -^ — 5-      , 


r2     2 
6-73 

Xy^ 


-?  +  P 


Y  =  -Y^  =  i 


c2      2        ' 

6-y3 


ti'   =  1  -  p 


2     5 


X^  =  i  (ati'-p')  +  i  J    (ati'+p')^  -  k[i'r 


ii)   Horizontal  Magnetic  Field:  y^  =  0 


5-X   » 
1/2 

t,8] 


(a)  Field  along  x-axis:  y^  =  Jo  =  0  ("north-south'  transmission 
in  magnetic  equatorial  regions): 
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{■6,9) 


X 
6 

X(5-X) 


a  =  l  -3 


p  =  1  - 


(e^-y^-xe) 


p '  =  -?  +  p' 


c  =  -c->  =  i 


pxy. 


(5^-yJ-X5) 


J  ti«    =  1  - 


\2  =  i  (cm,'-?')  ±  I  |(a^'+?')^-UacH 


2/^2     2^ 
P  (5  -y^) 

(6^-yJ-X5) 


(b)  Field  along  y-axis  =  y3_  =  y3  =  0   ("  east-west"    transmission  in 
magnetic  equatorial  regions)  ; 

r 


X(S-X) 


a  =  1  - 


(8.10) 


(e^-y^-xe) 


P   =  l-g 


a  =-a^i-  =  i 


J         pi   X.  _3  +  p' 


pxy. 


(e^-y^-xe) 


H'    =  1  - 


2/e2     2. 
p  (6  -y^) 

(e^-y^-xs)  ' 


r>  2 

X     =  a     +  ati',   and 

2 

B,     Normally  Incident  Plane  Waves;     p  =  q  =  0 
In  these  cases  A(z)  has  the  general  form 


(8.11)       a(z) 


1 
0 
0 
0 


with  the  characteristic  equation 

(8.12)  (X^  -  a)(^^  -  P)   -  rf'^  =  0» 

Thus  for  normal  incidence  the  characteristic  equation  is  always  biquadratic, 

and  the  eigenvalues  are  given  by 


(8.13) 


X^  =  \   (a+}3)  ± 


'{ 


(a-p)  +  UyY^ 


1/2 
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For  the  indicated  magnetic  field  directions  we  have: 

i)  Vertical  magnetic  field  =  y-,  =  yp  =  0     1.9],  [5]: 


(8.1U) 


Q   1    X6 
a  =  p  =  1  ---^ 

5-y3 

xy, 
Y  =  _Y^.  =  i  -^ 

6  -y. 


X  =  a  +  iyj 


ii)  Horizontal  magnetic  field  »  y^  =  0  ]9] ,    [5]: 

(a)  Field  along  x-axis  =  yp  "  y^  "  0: 


(8.15) 


Y  =  Y»-  =  0 
a  =  1  -  V6 

_  .            X(6-X) 
P  -  1 p     p 

X     =  a  and 
X^  =  Pi 

(b)  Field  along  y-axis  =  y,    "  y-3  =  0 


same  as  equations  (8.l5)  with  y^  replaced  by  y,    and  a 
and  p  interchanged. 


iii)  Oblique  magnetic  field  in  yz-plane  =  y,    =  0       [9] , [loj , [l3j : 

X5(6-X) 


a  =  1  - 


[(5-X)(52-y^)-6y^] 


.2     2     2 


(8.16) 


P  = 


[5  -y3-y2-X5](6-x) 
[(5-X)(52-y^)-6y^] 


Xy3(5-X) 

Y  =-yK-  =  i J5 — ^ ?5- 

[(6-X)(62-y^)-6yP 

X^  -  ^  (a+i3)   ±  ||(a-p)^-UY^ 


1/2 


-  29  - 


and 

iv)  Oblique  magnetic  field  in  xz-plane:  72=0* 

same  as  equations  (&.I6)  with  y^   replaced  by  y^ 
and  a  and  p  interchanged. 

9,  Normal  Incident  Wave  and  Oblique  Magnetic  Field 

The  solution  of  the  problem  formulated  as  case  B,  iii)  in  the  pre- 
ceding Section  shall  be  developed  here.  For  this  configuration  we  have 
p  =  q  =  0  (i.e.,  normal  incidence)  and  y^  =  0  (oblique  magnetic  field  in 
yz-plane).  For  the  longitudinal  field  components  we  have,  from  equations 
(8.2),  (2.2)  and  (2.lj): 


(9.1) 


13(2)  = 
V,(z)  = 


-Xy, 


[(5-X)(6^-y^)-5y^] 


[i5V^(z)+y3V2(z)]. 


The  transverse  components  must  satisfy 

1° 


('•^'     k 


.     CO 

1  — 
c 


(z)    / 


\ 


1 

0 
0 


0       0     \ 
0        Y 

0    -p 


0-10 


/ 


V^(z) 

I.;(Z) 

c 


where  the  elements  a(z),  p(z)  and  y(z)   ^^  given  explicitly  in  equations 
(8.16).     This  system  can  be  reduced  to  two  second-order  coupled  equations 
in  a  variety  of  ways  and  has  been  considered  from  that  viewpoint  by  various 
authors  [lL,13  ].     However,  we  use  here  the  method  of  Section  U. 

Thus  we  label  the  eigenvalues  of  the  coefficient  matrix,  given  by 
equations   (8.16),   as  follows: 


\ 


=  -\r 


={^./^)^ 


-r 


(9.3) 


=  -X,.  =^/^  .)/r^)W 


^3'  "\~Y    2        V'   2 
Let  us  introduce  the  quantity  u(z)  defined  as 
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.(.,  = .  (fe  .  Z^)^ 


(9.U) 


([ 


2y3(6-XJ 


t2 


2/3 (6-x; 


+  1 


Substituting  (9.U)  in  (9.3)  we  obtain  for  the  eigenvalues 


\- 


=  1  - 


(9.5) 


5+iy-u 


^3  '    5+iy3/u 


Since,   in  general,    [uj   ?<  1,  X^    /  0  and  \-,  ^  0,  the  eigenvalues  are  distinct 
and  the  coefficient  matrix  may  be  diagonalized  by  a  similarity  transformation 
(i.e.,  equation  (U.?)).     Such  a  transformation  is  determined  by 


/  W(zii 


P(z)  = 


1 

\ 
\ 

-1 


u(z)| 


(9.6) 


-u(z) 


P-^(z) 


2(l-u2) 


.-1  \ 


-1 


/ 


-U(Z); 


1 

-1 

s 

s 

\ 

s 

-1 

1 

^ 

-1 

^ 

1 

-t 

-1 

X-} 

1 

We  now  proceed  as  in  equations  ()4.9)  to  (U.I6)  to  find  the  solution. 
From  the  above  matrices  one  obtains 


-31- 


(9.7) 


'^i/^J 


p-1    dP 
'^        dz 


JlnJ)^^.l)    J^n^ 


(X^-Xj)"-^       -(\^+x^r^ 


-V  (d  in/^ 


s-1 


Ui  \  dz    J 


(x^-x^)'-"        (x^+x^) 


-1 


f/n/xju^-l)      /nyxl 


N-1 


(x^+x^)'-"        (x^-x^) 


-1 


i^nyx^  Injx^iu^^l). 


The  matrix  D(z)  of  eqioations  (4.11a)  has  as  diagonal  elements  the  four 
aiantities  d..,  where 


(9.8) 


d^^(z) 


d,2('!) 


^^-^(z) 


V^^) 


=     tl?V^)-^^nyx,(uSl 


=  t  i  ?  s  (^)  -  ^i-y  s(-'  - 1 


Thus^f^  )  D(i)    L  has  elements     exp 


d^^(x)  dx 


given  "by 
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(9.9) 


X,(z  )(u2(z  )  -1) 

1       O 0 


2     0 5 


ezp 


±i^/       X,(x)dx 


ezcp 


+  i  ^i  /     \,U)  dx 


—      c 


The  elements  of  the  niatrix  M(z)  =    Cm^.(z)    )  'become 


0         e 


12 


dz 


^21     ° 


xTcry 


V(z)  dinVn 


1^^^        [7x^(z)X3(.^ 


dz 


'11. 


X,-X3 


\    _!22 


(9.10) 


hk 


h^h 


i    _!2I.  _!2L. 


"^1       ^2 


"^21       "^22. 


^3       ^^ 


"'23      °'2^y 


X,(z^) 


V(z)  dins/u 


^1 


dz 


where  e.  .(z)  =     erp 


f. 


0         e 


34 


e,,^     0 


V'^3 


:^'[^iW-Xj(x)] 


«bc 


L22- 


x,-x  x^+x    I        .  \ 


^^1  ^^2 

^1*  ^^3  V  '^3/ 


Wi   °'42y 
W3    ""^y 


,  Using  the  definitions  (U.14) 


and  (4.15)  we  find  that  for  the  present  case  the  solution,  (^+.16),  takes  the  form 
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(9.11) 


X^(z^)(u  (z^)-l 
\^(z)(u^(z)-l 


e3p 


c^(z^,z) 
2 


u-(z) 


'X-(z^)(u  (z^)-l 

JO       0 

X^(z)(u^(z)-1 


exp 


+  i  ^/   X^U)  dx  1  c^(z^,z) 
^0        J   4 


This  solution  has  the  interpretation  given  in  Section  kB.  and  ^t.  The 

continuous  reflection  and  transmission  coefficients  S  (z),  defined  in  (4.1?) 

and  (4.18),  are  here  the  coefficients  of  the  e.  .  in  the  matrix  (9,10),  There 

are  essentially  four  of  these,  which  we  denote  by 


00 


-  d7  '^"V^i(z)  > 


(9.12) 


P^W)   - 


PQX       PlO 


^/n/1977, 


S^z;:  ^^n/d:7y 


yx^(z)x^(z) 


X^(z)  +  X^(z) 


T   =  t 
ox     xo 


vizl_  Ai,y7(7y 


X^(z)  -  X^(z) 


Here  p  (z)  and  p  (z)  are  the  ordinary  and  extraordinary  continuous  reflection 


00 


XX 


Coefficients  while  p  (z)  and  p  (s)  are  the  nixed  reflection  coefficients, 


ox 


xo 


The  "T:   (z)  and  T   (z)  are  the  continuous  transition  coefficients^ which  "by 


ox 


xo 


definition  are  of  mixed  type.  From  equations  (4.1?),  (9.IO)  and  (9*12)  we  have 


^2=^21=Poo,        \k  =  ^23  =  ^^(^o^/^l^'^o^   Pox. 


(9.13)     £3^,  =  E:,3  =  P^ ,        \^   =  R32  =  yy^TV^TTT   p^^  ^ 


^3  =  V  =  y^^V^Tx^:^  t:-,,; H3,  =  H,2  =  7^u,)/X3U,J    r,. 
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Thus  the  above  classification  is  seen  to  check  with  the  table  (U»l8), 

The  quantities  (9.12)  aid  in  simplifying  and  interpreting  the  continuous 

reflection  factors  c.(z  ,<z) , 

1  o' 

The  solution  (9 oil)  reduces  to  a  slight  modification  of  the  usual 

magneto-ionic  solutions  when  the  coupling  is  small  (or  neglected).  This 

result  is  obtained  by  assuming  the  c.(zq>e)  to  be  constants.  The  modifica- 

2     -1/2 
tion  consists  of  the  factor  (u  (z)-l)  '  ,  which  in  general  varies  to  the 

-1/2 
same  order  as  X  '  (i).  In  addition,  this  factor  clearly  shows  that  there 

2 

is  a  turning  point  (or  reflection  level)  -vAierever  u  (z)   =»  1, 

Polarizations,   defined  in  table  (U,22)   for  the   characteristic  waves 

are  evaluated  from  the  P(.z)  of  equation  (9o6):  are 

(9.1U) 


ORDINARY  WAVE 

-     -                            -         — , 
EXTRAORDINARY  WAVE 

Rising 

Downcoraing 

Rising 

Downcoming 

E-POLARIZATION 

-u(z) 

-u(z) 

-l/u(z) 

-l/u(z) 

H-POLARIZATION 

U(2) 

u(z) 

-l/u(z) 

-l/u(2) 

There  are  essentially  two  kinds  of  polarization,  the  ordinary,  u(z) 
and  the  extraordinary,  l/u(z).  These  quantities  were  first  derived  in  a 
more  direct  manner  by  Rydbeck  [.32] j  (see  his  work  for  a  more  detailed 
discussion) o  The  above  table  shows  that  the  electric  and  magnetic  polariza- 
tions are  130  out  of  phase.  However,  it  should  be  recalled  that  these  re- 
present the  polarizations  of  characteristic  waves  assumed  to  be  independently 
propagating  and  that  such  waves  are  not  necessarily  excited  by  a  source  out- 
side the  ionosphere, 

10,  Oblique  Incident  Wave  and  Vertical  Magnetic  Field 

In  exact  analogy  with  the  previous  section  we  shall  obtain  here  the 
solution  for  the  configuration  of  the  case  (A,i,)  in  section  8,  Thus  we  now 
have  y-i  =  yp  "  Q  "  0*  and  equations  (8,2)  become 
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(10.1) 


I^(z)  =  p  V^Cz) 

V  (2)  =  i  2.  Pl_   -J^ 


Here  we  have  eliminated  lAz)  "by  means  of  the  first  equation  of  the  following 
system,  which  determines  the  remaining  field  components: 


(10.2)     f 
dz 


/V^(z) 


l,U) 


\U) 


iw 
c 


\ 


0   |j.   0   0 


a   0   0   0 


V   0   0   P' 


^^iz) 


v^ 


l^iz) 


vyo, 


,0   0-1   0 


V.cz), 


\ 


Here  we  have  discarded  the  prime  on  the  p.  for  simplicity.  Equations  (8,3)  give 
the  elements  of  the  coefficient  matrix  e35)licite.  This  system  (10.2)  is  almost 
the  same  as  (9.2);  hov;ever,  the  sliglit  difference  will  he  seen  to  cause  a  much 
greater  complexity  in  the  present  results. 
The  eigenvalues  are 


h   =  -^2  =' 


u^-  g' 


4 


ua  •*•  P' 
2   y 


-  p.  V^ 


(10.3) 


X3  =  .^,=/  t^-^  V/^^)'-  .v^ 


We  introduce 


(lO.ij) 


u(^)  --  ^^^   */f^^)'- 


=  i<  - 


r    2„ 
p  73 

2(S-x) 


»■   2 

p  y- 


2(6-xT 


1  - 


2r  1 
P  6 


and  then  the  eigenvalues  hecome 

X(8-ly^u) 
1 


...-  ..3»,     2    2      X(6+iyu)       g 

(10.5)  xf  =  1  -  ^  -  P  ;  ^;  =  1  -  ^  -  P  TelocT 


(t^-vp 


i^'-vp 


1- 


6X 


(6^-y^) 


The  dlagonalizing  matrices 


u(z) 
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/'      1 


P(z)=! 


(10.6) 


I- 


and 


u(z) 


P-^.)  =        ^ 


2(^-u^) 


/       1 

-1 


X,         \ 


^l\  j 
jx/X 


Oi) 


u(z) 


(n) 


u(z) 


X^H 


3/ 


-1 


From  these  matrices  we  compute 


r,,(0        =     P-^  P 


-V  d  |?n  <AI 


/  (x^-x^r^       (x^+x^r-"-  \ 


=alf^  ^n  y^)| 
X^V    dz      J\ 


X_       dz 


(10.7) 


\^   (X^+X^)"^  (X^-X^) 


-1 


-u 


1         1        d  ^n  /g 

u    _^ 


^3     v^ 


1 

0 

1 
1 


(X^-Xj)"^       -(x^+x^)-^ 


(X,+X,)"^  (X,-X-)' 


'1  "3 


1  "3' 


/fin/x   (u%)         in  /xT 


l/nv^ 


in^l\^{M^-^)  I 


\i.M 


^    1 

0 

1 
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The  diagonal  elements  of  i  —  A   -  P~"TP  become 


(10.8) 


'22 


ii?\-^in/^^ 


ii2  xj.^in/q^?::;:^ 


33 

The  matrizant-^*-  |  D  A  of  this  diagonal  mstrix  has  the  elements 


•^  tt(z^)  /  Xi(2)[u2(z)  -  ti(z)]      r   °  ,  ^    J 


(10.9) 


X.(z)[u2(z)  -  n(z)] 


e^cp  [1  i  |/'x3(x)dx]  . 


From  equations  (10,7)  and  (10. 9)  the  matrix  M(z)  defined  by  (1;.13)  is 
obtained  as 

(10.10)        H(z)  -  ^\.   (z))  -  4-  ^h^     (S  )  , 

where 


"21 


-  P- 


V 


21 


'12 


z) 


/ 


13       Hi 


\  "®23 


'2U 
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3    -   -i-  /^ 
21  "  J^J   h^ 


®31    ®32 


®4l    ^kz 


22 


0    e 


«43   ° 


m^  »  m^  ,  m^2  -  ^2  > 


B>i2  "  >'^^  l^^^  "^2  ' 


•"21 


"^21  • 


The  m^   are  defined  as  in  (9,10)  and 


e.   = 


IJ 


ejqp 


f/      [X,(x)  -  X^(x)]  di 


The  solution  and  the  continuous  reflection   and  transition  coefficierts  may 
now  be  easily  obtained  as  in  Section  9.     It  should  first  be  noted  that  the 
expressions  of  the  present  Section  reduce  to  those  of  the  previous  Section  if 
we  set  n(2)=l«     This  was  to  be  expected,   since  the  systems  (9»2)   and  (10,2) 
then  become  formally  identical.     The  scattering  coefficients  are  all  sims  of 
essentially  the  previous  expressions  plus  terms  due  to  the  inhoraogeniety  of 
IJ-(z): 


oo 


(10,11) 


p     (z)  »  p     (z) 
ox  xo 


/^^    [,U,^^„/^_^^inv^J. 


From  the  P  matrix  of  equations  (10,6)  we  find  that  the  polarization 
table  (U,22)  now  becomes  explicitly 


(10.12) 
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nKDINAHY  VA7B 

EXTRAORDmARI  WAVE 

Rising 

Dovmcoaiing 

Hiaine 

Downcoming 

E-POLAEIZATION 

-u(z) 

-u(z) 

-p.(z)/u(z) 

-HJi(z)/u(z) 

R-POLASIZATION 

u(z)/li(z) 

u(z)/p.(z) 

l/u(z) 

l/u(z) 

Here  there  are  four  distinct  polarizations,  whose  relationship  to  each 
other  is  somewhat  surprising:  the  ordinary  E-wave  and  extraordinary  H-wave 
polarizations  are  the  negative  reciprocal  of  each  other,  and  similarly  for 
the  other  waves.  The   specific  dependence  of  there  polarizations  on  ionospheric 
properties  may  "be  ohtained  "by  using  equations  (10. i+)  for  u(z)  and  (8,8)  for 


-1*0  - 

APPENDIX  I 
Definitions  and  Identities  of  the  Modified  Peano-Baker  Method 

For  a  square  matrix  A(z)  with  bounded  integrable  elements 
a..(z),  the  matrizant   -TI   j  a(x) r  is  defined  by 


(I 


Z  Z  Z  X 

,1)    -^  [a(x)]  5  I  +  y  A(x)  dx  +  y  A(x2^  f       A(x^)dx^dX2  +  .< 


Z  Z            Z 

O  0            0 

Then  the  following  identities  can  be  proven: 

z               ^                   z  z 

(1.2)         -Q.    |d(x)|  =     exp[y  D(x)dxJ  »  (6^^  exp  [  /    d^(x)dx]    ), 


"z  "  "  z 

o  o 


where  D(x)  is  a  diagonal  matrix,  and 


z      I    z  T 


o  o''o(_o  o 

where  X  and  Y  are  square  matrices  of  the  same  order.     The  inverse 
of  any  matrizant  is 


z  T  z 

-1      <-  •^  o 


(I.li)         -0.       ^A(x)]    -     -A   j^A(z)/        . 

Another  useful  relationship  is 
z      -V       z 

(1.5)  -0  [^^^^f  =  -H-Ux))  -n.[A(x)/  . 
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